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Abstract 

A simple method is proposed to construct the spectral zeta functions required 
for calculating the electromagnetic vacuum energy with boundary conditions 
given on a sphere or on an infinite cylinder. When calculating the Casimir 
energy in this approach no exact divergencies appear and no renormalization 
is needed. The starting point of the consideration is the representation of the 
zeta functions in terms of contour integral, further the uniform asymptotic 
expansion of the Bessel function is essentially used. After the analytic con- 
tinuation, needed for calculating the Casimir energy, the zeta functions are 
presented as infinite series containing the Riemann zeta function with rapidly 
falling down terms. The spectral zeta functions are constructed exactly for 
a material ball and infinite cylinder placed in an uniform endless medium 
under the condition that the velocity of light does not change when crossing 
the interface. As a special case, perfectly conducting spherical and cylindrical 
shells are also considered in the same line. In this approach one succeeds, 
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specifically, in justifying, in mathematically rigorous way, the appearance of 
the contribution to the Casimir energy for cylinder which is proportional to 
ln(2vr). 
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I. INTRODUCTION 



A considerable achievement in theoretical investigations of the Casimir effect ]I],|2| was its 
calculation for massive fields (scalar and spinor) with boundary conditions on a sphere [|3],[|] . 
The various divergent contributions had been discussed in detail from the point of view of 
the general theory of adiabatic expansions (resp. heat kernel expansion). In a subsequent 
paper ]|] it was clarified in which cases the calculation of the Casimir energy, after the 
proper renormalization, yields a meaningful (unique) result and in which not independently 
of the regularization used. For a massive field a well defined result can be obtained in 
any case using the normalization condition proposed there. Instead, for a massless field 
the heat kernel coefficient 02 must vanish in order to allow for a meaningful calculation of 
the Casimir energy. For instance, this is the case for a material body characterized by a 
polarizability and a permittivity when the speeds of light inside and outside are the same or 
their difference is small. The vanishing of for the Dirichlet and Robin boundary conditions 
(and as a consequence for the conductor and bag boundary conditions) when taking inside 
and outside contributions together made Boyer's and all subsequent calculations possible 
and meaningful. When using a clever regularization (like the zetafunctional one ^W) it 
is even possible to avoid the appearance of divergencies other than that in the Minkowski 
space contribution at all. 

Practically every problem on calculation of the Casimir energy (or force) has been con- 
sidered multiply with employment of more and more effective and elaborated mathematical 
methods. For example, the first calculation of the Casimir energy of a perfectly conducting 
spherical shell carried out by T.H. Boyer in 1968 || has required computer calculations dur- 
ing 3 years |J. Later this problem was considered in many papers |10-12|. By making use of 



the modern methods JT3[ it can be solved practically without numerical calculations (with 
a precision of a few percent). It requires only the application of the uniform asymptotic 
expansion for the Bessel functions. 

In recent papers [|l4| , |T5f the Casimir energy of a compact ball |16[ and infinite cylinder 
has been calculated by making use of the mode-by-mode summation technique. In these 
problems two sums appear, over the roots of radial frequency equation at fixed value of 
angular momentum and then over angular momentum. The either of these sums is divergent. 
In papers [ |T^| , P~5[ ] for each of these summation a separate regularization has been used. 
The first summation was carried out by applying the contour integration with subsequent 
subtraction of the contribution of an infinite homogeneous space. The second sum was 
evaluated by making use of the Riemann ( function technique. However the procedure of 
analytic continuation, required by rigorous approach, has not been considered there. 

The present paper pursues the aim to eliminate the minor points of preceding considera- 
tions, i.e., the Casimir energy for two configurations mentioned above will be calculated by 
the rigorous ( function techniques, and the analytic continuation of the relevant spectral ( 
functions will be carried out exactly. An essential advantage of this regularization procedure 
is that no manifestly divergent expressions arise in its framework, and it gives final finite 
result without any subtractions (renormalizations). 

The layout of this paper is as follows. In Sec. II, the spectral zeta function is constructed 
for a compact ball placed in uniform endless medium when the light velocity is the same 
inside and outside the ball. As a special case the zeta function for perfectly conducting 
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spherical shell is also considered. In Sect. Ill the spectral functions for infinite cylinder are 
constructed under the same conditions. These results provide a firm footing for the previous 
calculations of the Casimir energy for given boundary conditions by making use of a "naive" 
zeta function method. In Sect. IV the obtained results are shortly discussed. 



II. CASIMIR ENERGY OF A COMPACT BALL UNDER THE CONDITION 

£// = 1 

In the £ function method ]7||| the Casimir energy Eq is defined in the following way. 
Let Wp's be the eigenfrequencies of the quantum field system under the influence of the 
boundary conditions, and let tD p 's be the same frequencies when the boundaries are removed. 
By making use of this spectrum one defines the ( function for the problem in hand 

C(s)=£«-S s )- (2-1) 

W 

Here the summation (or integration) should be done over all the quantum numbers {p} 
specifying the spectrum. The parameter s is considered at first to belong to region of the 



complex plane s where the sum ( |2.1| ) exists. Further the analytic continuation of ( |2.1|) to 



the point s — — 1 should be constructed. After that one puts 

E c = \as = -l). (2.2) 

Let us consider a solid ball of radius a, consisting of a material which is characterized by 
permittivity S\ and permeability fii. The ball is assumed to be placed in an infinite medium 
with permittivity e 2 and permeability /i 2 . The eigenfrequencies of electromagnetic field for 
this configuration are determined by the frequency equation for the TE-modes |T7] 



Aj (au) = yjE\Hi s'i(kia) ei{k 2 a) - ^e 2 Hi si(kia) ej(/c 2 a) = 0, (2.3) 
and the analogous equation for the TM-modes 



TM, 



auj) = y/eopi s'^ha) ei(k 2 a) - ^/e\[i 2 si(kia) e[(k 2 a) = , (2.4) 



where ki = ^JSiHi u, i = 1,2 are the wave numbers inside and outside the ball, respectively. 
Here 5/(x) and ei(x) are the Riccati-Bessel functions 



vr , . s _ . s it 



si{x) = \J^: J i+i/2(x) , e,(ar) = yj ^ H^ 1/2 (x) , (2.5) 
and prime stands for the differentiation with respect to their arguments, k±a or k 2 a. The 



orbital momentum I in Eqs. (|2.3| ) and (2^) assumes the values 1, 2, . . .. 



As usual when one is dealing with an analytic continuation, it is convenient to represent 



the sum (|2.1| ) in terms of the contour integral 



>/n 2+1 / i/2 2N , /2 d , Aj E (az)A™(az) , s 

Cc(s) = —— lim & dz{z 2 + /i 2 " 7 V ln a tp a tm , 2.6 
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where the contour C surrounds, counterclockwise, the roots of the frequency equations in 



the right half-plane. For brevity we write in (2J3) simply A;(oo) instead of lim^oo Ai(az 



Transition to the complex frequencies z in Eq. ( |2.6| ) is accomplished by introducing the 
unphysical photon mass // 

u^{z 2 + ^) s l\^. (2.7) 

Extension to the complex z-plane of the frequency equations Aj E (az) and Aj M (az) should 
be done in the following way. In the upper (lower) half-plane the Hankel functions of the first 
(second) kind H^(az) (HjP(az)) must be used |18|] . Location of the roots of Eqs. ( |2.3| ) and 
( p.4|) enables one to deform the contour C into a segment of the imaginary axis (— iA, zA) and 
a semicircle of radius A in the right half-plane. When A tends to infinity the contribution 
along the semicircle into (baii{s) vanishes because the argument of the logarithmic function 
in the integrand tends in this case to 1. As a result we obtain 

„ , , (2i + l)o«,. /+*» , 2 , 2W2<<, <V E (z)A™(z) 

Now we impose the condition that the velocity of light inside and outside the ball is the 

same 

£i/ii = e 2 ^2 = c~ 2 . (2.9) 
Under this assumption the argument of the logarithm in ( |2.8| ) can be simplified considerably 



[ Hfl with the result 



\aj r^i n Jo ay 



where 



^J 1 -^, °i{y) = ^-[si{y)ei{y)\. (2.11) 
E\ -+- fc2 cty 

Here si(y) and e;(y) are the modified Riccati-Bessel functions 




si(x) = \j—Iu{x), e l {x) = \\—K v {x), v = 1 + 1/2. (2.12) 

More details concerning the contour integral representation of the spectral ( functions can 
be found in [|,|ra|-|5T|] . 

Further the analytic continuation of Eq. ( |2.10| ) is accomplished by expressing the sum 



over / in terms of the Riemann ( function. This cannot be done in a closed form. Making 
use of the uniform asymptotic expansion (UAE) for the Bessel functions in increase powers 
of \ jv enables one to construct the analytic continuation looked for in the form of the series, 
the terms of which are expressed through the Riemann ( function. The problem of the 
convergence of this series does not arise because its terms go down very fast. 
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We demonstrate this keeping only two terms in UAE for the product of the Bessel 
functions I u {vz)K,.(vz) l2l 



IJvz)KJvz) ~ 



t 

2^ 



t 2 (l - 6t 2 + t 4 ) 

1 + — - + ••■ 



1 



(2.13) 



After changing the integration variable y — vz in Eq. ( |2.10|) we substitute Q2.13| ) into this 
formula and expand the logarithm function up to the order u~ 3 keeping at the same time only 
the terms linear in £ 2 . The last assumption is not principal. It is introduced for simplicity 
and in order to have possibility of a direct comparison with the results of Ref. ||14|1 . Thus 
we have 




1 2 



sr 8 r s 

2 v 2 16z/ 4 



^-(zI v (vz)K v (vz)) 
dz 



■12 + 216t 2 - 600£ 4 + 420t 6 ) + 0(v~ 



(2.14) 



Integration over z can be done by making use of the formula 



1 
2 



r - 



[2.15) 



Also the properties of the T function 
T(z)T(l-z) 



7T 



T{l + z) = zT{z) 

sirnrz 

prove to be useful. After simple calculations we arrive at the result 

e 



(ball(s) 



32 \a 



s(2 + s )(4 + s) [y: v- i - s + p(s) e »~ 3 ~ s + 



(2.16) 



(2.17) 



i=i 



v 



1 + 1/2 



where 



p(s) 



l-- 4 (6 + s) + ^6 + s)(8 + s) 



192 



(6 + s)(8 + s)(10 + s) 



(2.18) 



The zeta function Cbaii( s ) represented in the form (|2.17|) is defined for Re s > due to the 
first sum over I. This term corresponds to the order 1/v in the uniform asymptotic expansion 
( 2.13 ). The second sum in ( |2.17| ), defined at Re s > —2, has been generated by the term 



1/z/ 3 in Eq. ( [2.13] ). It is clear that the terms of order l/v 2k+1 in (|2.13|) will give rise to the 
singularity of Cbau{ s ) a t the points s = —2k, k = 0, 1, 2, . . .. Due to the multipliers in front 



of the square brackets in (|2.17|) the first three singularities are really the indefinitenesses like 
• oo. 
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The analytic continuation of Eqs. Q2.17|) , (|2.18|) into the region Re s < can be accom- 
plished by expressing the sums over angular momentum I trough the Riemann ( function 



according to the formula |23 



1=1 



{T - l)C(s) 



v = 1 + 1/2. 



(2.19) 



As a result one gets 



(ball(s) 



e (o 



s{2 + s)(4 + s){(2 1+s - 1)C(1 + s) - 2 



)1 + S 



^2.20) 



32 \o. 

+^)[(2 3+s -l)C(l + ^)-2 3+ 1 + ...}. 

The singularities in Eq. (|2.17j) are transformed in ( 2.20|) into the poles of the Riemann ( 
functions at the points s = 2k, k = 0, 1, 2, . . . 

1 



C(l + s) ^ ~ + 7 + 
C(3 + s) ^ 7^ + 7 + --- 
C(5 + s) ~ j^ + 7 + --- 



-2 

-4. 



(2.21) 



where 7 is the Euler constant. The first three poles are annihilated by the multipliers in 
front of the curly brackets in Eq. ( J2.2U ) . The first surviving singularity (simple pole) appears 
only at the point s = —6. Thus the formula Q2.20 ) affords the required analytic continuation 
of the function (baii(s) into the region Re s < 0. In view of Eq. ( |2.2j ) we are interested in 
the point s = — 1 where Cbaii(s) given by ( |2.20| ) is regular 



(ball( — 1) 



3^c 
32a 



9 7T 2 \ 

1 + 4 +... 

128 V 2 / 



(2.22) 



It is exactly the first two terms in Eq. (3.10) of Ref. [|14 |. The procedure of analytic 
continuation presented above can be extended in a straightforward way to the arbitrary order 
of the uniform asymptotic expansion (|2.13|) . Certainly in this case analytical calculations 
should be done by making use of Mathematica or Maple. 

The problem under consideration with £ = 1 is of a special interest because in this case 
it gives the Casimir energy of a perfectly conducting spherical shell. As it was noted above, 
this configuration has been considered by many authors. We present here the basic formulae 
which afford the analytical continuation of the corresponding spectral ( function. We again 
content ourselves two terms in the UAE (|2.13|) . In the next formula Q2.14 ) It is impossible 
to put simply £ = 1 in the next formula Q2.14D - One has to do the expansion here anew 
keeping all the terms ~ 1/zA This gives 




az 



(2.23) 
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_L^ 8 + 3 zf8 r_ x + lgf2 _ 5()t4 + 35t6 x + e 

After integration and elementary simplifications we arrive at the following result for the 
spectral function in hand 



(shell(s) 



where 



32a- 



q(s) 



s(2 + s)(4 + s) 



E- 



-1-s 



+ q{s) E ^ 3 " S + • • • 



l=i 



3840 



(480 + 868s + 504s 2 + 71s 3 



(2.24) 



(2.25) 



Obviously formula (|2.23j ) has the same singularities as Eq. ( [2.17| ), i.e., it is defined for 
Re s > 0. The analytic continuation is accomplished by making use of Eq. ( [2.19| ) 



(shell{s) 



32a- 



s(2 + s)(4 + s){(2 



l+s 



l)C(l + s)-2 1+s + 



(2.26) 



+g(s)[(2 3+s -l)C(3 + s)-2 3+s ] + ...} 



The nearest singularity in this formula is simple pole at s = —6. As above it is originated in 
the term ~ 1/u 7 in the UAE ( p. 13 ). At the point s = —1 the spectral zeta function ( s heii{s) 
is regular and gives the following value for the Casimir energy of a perfectly conducting 
spherical shell 



E 



shell 



-i) = ttC 



shell \ 



-1) 



64a 



7T 



256 



-4 +. 



0.046361 



(2.27) 



Without considering the analytic continuation and do not carrying out the analysis of the 



singularities in the complex s plane this result has been obtained in [13 



Undoubtedly, the calculation of the Casimir energy of a nonmagnetic dielectric ball 
(eifii 7^ £2/^2) by a rigorous ( function method is also of a special interest. However, in this 
case the very definition of the spectral zeta function should probably be changed in order 
to incorporate the contact terms which seem to be essential in this problem 



III. VACUUM ENERGY OF ELECTROMAGNETIC FIELD WITH BOUNDARY 
CONDITIONS GIVEN ON AN INFINITE CYLINDER 



Calculation of the Casimir energy of an infinite cylinder p8| , p!5| proves to be a more 



involved problem in comparison with that for sphere. In this section the spectral zeta 
function Ccyi{s), for this configuration will be constructed its analytical continuation into 
the left half-plane of the complex variable s will be carried out, and relevant singularities 
will be analyzed. 

Thus we are considering an infinite cylinder of radius a which is placed in an uniform 
unbounded medium. The permittivity and the permeability of the material making up 
the cylinder are S\ and fix, respectively, and those for surrounding medium are e 2 and 
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H2- We assume again that the condition ( |2.9|) is fulfilled. In this case the electromagnetic 
oscillations can again be divided into the transverse-electric (TE) modes and transverse- 
magnetic (TM) modes. In terms of the cylindrical coordinates (r, 9, z) the eigenf unctions of 
the given boundary value problem contain the multiplier 

exp (±iut + ik z z + in9) (3.1) 

and their dependence on r is described by the Bessel functions J n for r < a and by the Hankel 
functions or for r > a. The frequencies of TE- and TM-modes are determined, 
respectively, by the equations |lj 



A™(A, a) = \a[inJ' n {\a)H n {\a) - ^ x J n (Xa)H n (Xa)'\ = , (3.2) 
A™(A, a) = \a[sM\a)H n (\a) - e x J n (Xa)H n (Xa)'} = , (3.3) 

where n = 0,±1,±2, .... Here A is the eigenvalue of the corresponding transverse 
(membrane-like) boundary value problem 

X 2 = ^~kl (3.4) 

In a complete analogy with the preceding Section we define the Casimir energy per unit 
length of the cylinder trough the spectral zeta function 

E cyl = ^Ccyi(-l) ■ (3.5) 

Let X nr be the roots of the frequency equations (|3.2|) and (|3.3| ), then the function ( cy i(s) is 
introduced in the following way 

/+oo rjh 
-± £[(A nr (a) + k 2 z )- s !* - (A nr (oo) + k 2 )^ 2 } . (3.6) 
-oo Z7T n r 

In terms of the contour integral it can be represented in the form 

Again we can take the contour C to consist of the imaginary axis (+ioo, —ioo) closed by a 
semicircle of an infinitely large radius in the right half-plane. Continuation of the expressions 
A™(Aa) and A™(Aa) into the complex plane A should be done in the same way as in the 
preceding Section, i.e., by using H^(X) for Im A < and H^(X) for Im A > 0. On the 
semicircle the argument of the logarithm in Eq. ( |3.7|) tends to 1. As a result this part of the 
contour C does not give any contribution into the zeta function ( C yi{s)- When integrating 
along the imaginary axis we choose the branch line of the function 0(A) = (A 2 + k 2 z )~ s l 2 
to run between —ik z and +ik z , where k z = > and use further that branch of this 

function which assumes real values when \y\ < k z , where y = Im A. More precisely we have 

e -W2( y 2 _ k 2 z )~ s / 2 , y > k z , 
<j>(iy) = { (%-y 2 )-' 2 , \y\<k z , (3.8) 

2 ™/2(y2 - k 2 z )~ s/2 , y < -k g . 
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Employment of the Hankel functions H^(X) and H^(X) by extending the expressions 
A^ E (A) and A™ (A) into the complex plane A, as it was noted above, gives rise to the 
argument of the logarithm function depending only on y 2 on the imaginary axis. It means 
that the derivative of the logarithm is odd function of y. As a result the segment of the 
imaginary axis (—ik z , +ik z ) gives zero, and Eq. (|3.7Q acquires the form 



(cyl(s) 



C ITS 
7T Z I 



+°° poo poo 

£ dkj (y 2 - k 2 z y s / 2 d y \n 

,-_^J0 Jk z 



A™(gy) n A™(ay) 

A™(ioo)A™(lOo) • 



(3.9) 



Changing the order of integration of k z and y and taking into account the value of the 
integral 



jyuv'-kir'^^y 1 



r 1 



Re s < 2 



(3.10) 



we obtain after the substitution ay — > y 

1 



c 



cyl\ 



-s +oo 



E 



d 



dyy^--\n[l-e^ n (y)} 
u dy 



(3.11) 



where 



(*n(y) = y{ I n{y)K n (y))\ £ 



£i + ^2 



(3.12) 



We shall again content ourselves with the first two terms in the uniform asymptotic expansion 
( ^.13j ) and take into account only the terms linear in £ 2 . In this approximation, upon 
changing the integration variable y = nz, n = ±1, ±2, . . ., we have 



in ii-e 



z^-(I n (nz)K n {nz)) 
dz 



(3.13) 



zH 6 

An 2 



An 2 



(3 - 30i 2 + 35t 4 ) + 0(n 



Now we substitute ( |3.13| ) into all the terms in fl3.11| ) with n ^ 0. The term with n = in 
this sum will be treated by subtracting and adding to the logarithmic function the quantity 



e 



4 (l + y 



2\3 



(3.14) 



As a result the zeta function ( C yi(s) can be presented now as the sum of three terms 



Cci(s) = Z^s) + Z 2 (s) + Z 3 {s) 



(3.15) 



where 
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ZAs) 



Zo(s) 



2^aT ( - ) T 



-r - 



3-s 



dyy 



l-s 



d f 

dy \ 



2j2n- s - 1 + 1 



n=l 



3-s 



ln[l + 



d 



Z 3 (s) 



-s +00 

En" 3 - 

n=l 



32v^Far(-)r 



3-s 



^ z i-^[/t8( 3 _30t 2 + 35t 4 )] 



dz 



(3.16) 



(3.17) 



(3.18) 



In these equations Z\(s) has accumulated the term with n = from Eq. ( |3.11[ ) subtracted 
by (p. 14 ); Z 2 (s) involves the contribution of the term of order 1/n 2 in expansion (|3.13| ) and 
the added expression ( |3.14| ); Z$(s) is generated by the terms of order 1/n 4 in the expansion 

Taking into account that 

1 3 

Mo(2/)l^o^l and #j(3/)ly-oo -> ^ + ^-j , (3.19) 
the integration by parts in Eq. (|3.16 ) can be done for —3 < Re s < 1 with the result 



ZAs) 



s- 1 



2vW(-)r 



3-s 



ay Jo 



dyy~ s L[l-e Hl(y)} + jy 4 t 6 (y)\ . (3.20) 



With allowance for ( |3.19| ) one infer easily that the function Z\(s) is an analytic function of 
the complex variable s in the region —3 < Re s < 1. In the linear order of £ 2 it reduces to 



z[ m (s) = r 



s - 1 



2^aT ( - ) T 



3-s 



dyy~ 



4(1 + y 2 ) 



(3.21) 



This function is also analytic in the region —3 < Re s < 1. Integration in Eq. ( |3.17| ) can be 
accomplished exactly by making use of the formula 



r°° d 
This gives for Z 2 (s) in ( |3.17| ) 



s-l r /l + s 



5- s 



Z 2 (s) = £ 



2 , c\- s (1 - s)(3 - s) 



64y / 7r a 



2E 



7? 



-s-1 



+ 1 



n=l 



1 < Re s < 5 . 



1 + s 



(3.22) 



(3.23) 



In view of the sum over n in ( |3.23| ) the function Z 2 (s) is defined only for s > 0. 
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For simplicity we apply in Eq. (|3.18 ) the integration by parts which is correct for —3 < 
Re s < 2 and leads to the result 



z 3 (s) = e 



c\- s (1 - s)(3 - s)(7s 2 - 4s - 27) 
a, 



3 + s 



61440ra 



n=l 



n 



-s-3 



(3.24) 



Again the sum over n in ( |3.24j ) gives the restriction Re s > —2 for definition of the function 
Z 3 (s). 

Thus the spectral zeta function ( C yi{s) in the linear approximation with respect to £ 2 and 
with allowance for the first two terms in the UAE (|3.13|) is given by 



C cy i(s) = Zf n (s) + Z 2 (s) + Z 3 (s) 



(3.25) 



where the Z's are presented in Eqs. ( |3.21| ), ( |3.23| ) and Q3.24| ), respectively. Summing up 
all the restrictions on the complex variable s which have been imposed when deriving Eqs. 
( |3.21| ), ( |3.23| ), and (|3.24|) , we infer that ( cy i(s) is defined in the strip < Re s < 1. In order 
to continue these equations into the surroundings of the point s = —1, it is sufficient to 
express the sum in Eq. (|3.24j ) in terms of the Riemann ( function 



z 2 (s) = e 



c\' s (l-s)(3-s) 
64a/7t a 



r 



[2C(s + l) + l]- 



l+s 



(3.26) 



It is left now to take the limit s — > — 1 in Eqs. ( |3.21| ), ( ^.23| ) and ( ^.26| ). A special care should 
be paid when calculating this limit in ( |3.26| ) in view of the poles of the function T((l + s)/2) 
at this point. Using the values 



C(o) 



l 



, C'(o) 



l 



in2vr, r'rr 



7: 



one derives 



lim [2C(1 + s) + l]r 



s-»-l 



1 + S 



= hm [2C(0) + 2C'(0)(1 + s) + 0((1 + s) 2 ) + 1] 

s— >— 1 

With allowance for this we obtain from ( |3.26| ) 

ce i 



+ 7 + 0(l + s) 



Z, -1 



2vra 2 4 



-ln(27r) . 



(3.27) 

(3.28) 
-21n(27r) . 

(3.29) 



The appearance of the finite term proportional to ln(27r) is remarkable for the problem 
under consideration. It is derived here in a consistent way by making use of an analytic 
continuation of the relevant spectral zeta function. In Ref. [|15[ it was obtained in a more 
transparent though not rigorous way. 

Gathering together Eqs. ( ET2TD , ( |3T24|) with s = -1 and Eq. ( |3729| ) we have 
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(cyl(-l) 



ce 

2na 2 

ce 

2na 2 

ce 

2na 2 



yd y 



V 1 



4(1 + y 2 ) 3 
(-0.490878 + 0.034269 + 0.459469) 

0.002860. 



1 +°° 1 1 
+ — T 4t + - ln(27r) 
48 ^n 2 4 v ' 



(3.30) 



This result is not the final answer in the problem in hand. The point is that in view of 
severe cancellations in ( |3.30| ) the contribution of the next term in the UAE ( |3.13| ) proves to 
be essential. Its account gives 



Ccyl(-l) 







(3.31) 



Thus the Casimir energy of a compact cylinder possessing the same speed of light inside 
and outside proves to be zero. The consideration presented in this Section can be extended 
to the next term of order ~ 1/n 6 in the UAE ( |3.13| ) in a straightforward way. Therefore we 
shall not present here these rather cumbersome expressions [ j27j] . 

Now we address to the consideration of a special case when £ = 1. It corresponds to a 
perfectly conducting cylindrical shell | 15| . Instead of the expansion (|3.13| ) we have 



z^-(I n (nz)K n (nz)) 
dz 



(3.32) 



4 + 6 



zH 



An 2 



1 + ^ ( 3 - 30 ^ + 35 * 4 + \ zH ') + °(™~ 4 ) 



Proceeding in the same way as above we obtain for the spectral zeta function concerned 



*cyl 



(3.33) 



where Zi(s) is given by Eq. ( |3.20| ) with £ = 1, Z 2 (s) is the same as in Eq. (|3.26| ), and Z 3 (s) 
now is 

'3 + 



ZJs] 



;i - s)(3 - s)(71s 2 - 52s - 235) 



r 



61440v/7ra 



l-s 



E 

n=l 



n 



-3-s 



(3.34) 



At the point s 



*shell 
scyl 



1 it has the value 
1 



-r 



2na 2 
1 

1 



0.6517) + 



2<ko? 480 ^ n? 



+ 



8vra 2 



ln(27r) 



-0.6517 + 0.0240 + 0.4595) 



0.0268 



(3.35) 



This exactly reproduce the contribution of the first two terms in calculations of the Casimir 
energy for cylindrical shell in Ref. [TH] . With higher accuracy this energy is given by 
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E^f = — -0.01356. (3.36) 



In a recent paper ]29] the vacuum energy of a perfectly conducting cylindrical surface 
has been calculated to much higher accuracy by making use of another version of the zeta 
function technique. By integrating over dk z directly in Eq. ( |3.6| ) the authors reduced this 
problem to investigation of the zeta function for circle, which has been considered earlier 
by introducing the partial wave zeta functions for interior and exterior region separately. In 
this respect our approach dealing only with one spectral zeta function for given boundary 
conditions proves to be more simple and straightforward. 



IV. CONCLUSION 

The method for constructing the spectral zeta functions developed here proceeds from 
the contour integral representation with a subsequent employment of the uniform asymptotic 
expansions for the Bessel functions. Upon an analytic continuation the zeta functions prove 
to be presented as (infinite) series over the Riemann ( functions with rapidly decreasing 
terms (see, for example, Eqs. ( |2.20|) , ( [2.26|) ). 



We didnot pursue here the goal of obtaining high accuracy when calculating the Casimir 
energy. In fact we seek to present the consideration in such a form that no manifest diver- 
gencies appear. An obvious advantage of the regularization method in hand does not need 
any renormalization. 

By treating the boundary condition given on an infinite cylinder, we have clearly demon- 
strated the importance of a consistent analytic continuation of the relevant spectral zeta 
function, in contrast to identifying simply the sum of the type Yl^Li n ~ s with Riemann ( 
function, in order to involve correctly the contributions to the Casimir energy proportional 
to ln(27r). 

Consideration in this framework of the same configuration of vacuum electromagnetic 
field but with different velocities of light inside and outside the boundaries probably will 
demand the modification of the definition of the spectral zeta functions for incorporating in 
a proper way the contact terms important in this case p4|-26 . 
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